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摘要 : 采用 2019 年 10 月 至 2020 年 3 月 北斗 观测 数据 获得 单 天 台 站 坐标 ， 实 现 BDCS 并 与 
ITRF2014 对 齐 。 利 用 单 GPS 观测 数据 的 解 算 台 站 坐标 和 IGS 提供 的 台 站 坐标 对 BDS 观测 数 
据 解 算 结 果 进 行 精 度 评估 ， 比 较 结 果 表 明 ， 采 用 BDS 数据 解 算 的 台 站 坐标 与 IGS 提供 单 天 台 站 
坐标 差异 的 RMS 均值 小 于 1 cm， 与 单 GPS 数据 解 算 台 站 坐标 差异 的 RMS 均值 小 于 1.5 cm. 
BDS 数据 解 算 台 站 坐标 时 间 序 列 弥散 度 偏 大 ， 并 对 其 原因 进行 了 初步 的 分 析 。 同 时 采用 GPS 
测 数 据 初 步 估计 中 国 区 域 监测 站 速度 场 ， 随 着 日 后 北斗 导航 系统 的 稳定 运行 以 及 北斗 观测 数 
据 的 累积 ， 有 望 获得 更 高 精度 的 台 站 坐标 和 速度 及 BDCS 向 ITRF 的 对 齐 。 
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1 3 言 


2020 年 7 月 31 日 ， 北 斗 卫 星 导 航 系 统 (BeiDou Navigation Satellite System, BDS) FF 
始 向 全 球 及 中 国 周边 地 区 提供 7 种 类 型 服务 "， 其 中 通过 卫星 播发 的 导航 服务 包括 : 向 全 球 
提供 米 级 的 导航 定位 授时 服务 ， 向 中 国 及 周边 地 区 提供 的 分 米 级 星 基 增强 服务 和 分 米 至 厘 
米 级 精密 单 点 定位 服务 ””。 与 其 他 GNSS 系统 不 同 ， 北 斗 三 号 地 面 监控 站 主要 分 布 于 中 国 
境内 ， 监 控 站 地 理 分 布 无 法 覆盖 全 球 ， 会 对 卫星 定 轨 、 授 时 和 全 球 服务 产 生 影响 ， 因 此 ， 星 
上 搭载 星 间 链 路 设备 ， 利 用 星 间 双向 测量 数据 分 离 卫 星相 对 钟 差 和 相对 几何 距离 解 耦 卫星 
轨道 和 钟 差 ， 再 把 星 间距 离 作为 观测 量 结合 地 面 测量 数据 进行 星 地 星 间 联合 定 轨 ， 由 此 实 
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现 境 外 测量 “。Chen 等 人 ”对 北斗 三 号 基本 导航 服务 性 能 进行 评估 的 结果 表明 ， 在 有 星 
间 链 路 数据 情况 下 ， 广 播 星 历 用户 距 离 误差 优 于 0.1 m， 广 播 钟 差 参 数 预 报 精 度 优 于 1.5 ns, 

空间 信和 号 精度 优 于 0.6 m 
导航 系统 的 高 精度 稳定 服务 以 及 兼容 互 操作 性 能 都 依赖 于 导航 系统 参考 框架 实现 和 维 
持 精度 。 通 常 卫 星 导 航 系统 采用 地 面 监测 站 接收 的 观测 数据 进行 卫星 精密 定 轨 与 时 间 同 步 
处 理 ， 生 成 导航 电文 信息 ， 实 现 全 球 导 航 定位 授时 服务 。 而 导航 系统 空间 基准 建立 与 维持 的 
目的 是 精确 测定 并 预报 地 面 监测 站 在 统一 的 地 球 框架 中 的 位 置 和 速度 ， 有 具体 方法 是 利用 空 
间 大 地 测量 方法 获得 监测 站 历 元 坐标 精确 估计 ， 再 利用 时 间 序 列 分 析 方 法 对 监测 站 运动 进 
行 建 模 ， 估 计 监 测 站 运动 的 线性 和 非 线 性 模型 参数 ， 实 现 地 球 参考 框架 的 维持 。 监 测 站 历 元 
坐标 估计 精度 与 卫星 轨道 精度 和 测量 模型 有 关 ， 参 考 框架 维持 精度 与 监测 站 运动 建 模 精度 
I Ke 
目前 ， 四 大 卫星 导航 系统 (GPS, GLONASS, Galileo 和 BDS) 均 采 用 卫星 大 地 测量 技 
术 实 现 监 测 站 坐标 估计 ， 采 取 定 期 向 国际 地 球 参 考 框架 (International Terrestrial Reference 
Frame, ITRF) 对齐 的 方法 实现 空间 基准 的 高 精度 和 系统 间 兼 容 互 操作 性 。 由 于 各 导航 系 
统 的 地 面 站 分 布 范围 、 建 站 和 条件、 观测 方法 等 各 不 同 ， 导 航 系 统 空间 基准 实现 和 维持 的 具体 
方法 存在 一 定 差异 。 美 国 GPS 坐标 系 WGS84 得 益 于 台 站 全 球 均 匀 分 布 和 GPS 卫星 成 熟 
的 测量 模型 及 动力 学 模型 。 全 球 均匀 布 站 易于 消除 区 域 性 整体 运动 对 全 球 参考 框架 维持 精 
度 的 影响 ， 成 熟 的 GPS 处 理 方法 包括 统一 系统 参数 、 精 化 天 线 相位 中 心 模型 等 ， 完 备 的 参 
考 框架 处 理 方法 最 终 实现 WGS84 与 ITRF 的 高 精度 对 齐 和 长 期 维持 ，WGS84 (G1762) 与 
ITRF2008 的 差异 为 1 cm ” ”。 Galileo 坐标 系 GTRF20V02 (Galileo Terrestrial Reference 
Frame) 与 ITRF2014 对 齐 精度 小 于 1cm ” ”。GLONASS 坐标 系 PZ-90 (Parametry Zemli 
1990) 对 齐 ITRF2014 精度 小 于 2 cm 
从 2013 年 开始 ， 北 斗 卫 星 导航 系统 采用 2000 中 国 大 地 坐标 系 (CGCS2000)”“”。 随 着 
北斗 全 球 系统 的 建设 ， 北 斗 导航 系统 服务 范围 从 区 域 扩展 到 全 球 ， 从 2017 年 12 月 开始 ， 北 
斗 导 航 系统 开始 采用 北斗 坐标 系 (BeiDou Coordinate System, BDOS)”. TIEA A 
Wu” 对 北斗 监测 站 分 别 进行 了 4 期 单独 的 GPS 观测， 利用 64 个 全 球 均匀 分 布 IGS 台 站 和 
8 个 中 国 区 域 台 站 的 GPS 观测 数据 ， 首 次 实现 BDCS 并 对 齐 ITRF2014， 监 测 站 坐标 精度 
= 为 1 cm， 速度 为 2 mm/a. Hi A?" FA IGS/MGEX 全 球 观测 网 台 站 GPS 与 BDS 数 
- TK GET SN, FR REAG AA ER TIE, KURA Ew, KAÉ 
BDS 数据 是 北斗 二 号 观测 数据 。 北 斗 二 号 系统 主要 为 亚太 地 区 用 户 提供 定位 、 测 速 、 授 时 
和 短 报 文通 信服 务 ， 在 轨 卫 星 数 为 5GEO/7IGSO/3MEO， 而 且 主 要 履 盖 区 域 仅 亚太 地 区 。 
随 着 北斗 三 号 全 球 服务 开通 ， 在 轨 卫 星 已 增加 至 3GEO/3IGSO/24MEO， 卫 星 履 盖 区 域 也 
扩展 至 全 球 ， 因 此 利用 北斗 单 系统 卫星 数据 实现 北斗 坐标 系 的 可 行 性 和 精度 有 待 进一步 分 

析 研 究 。 

2019 年 12 月 ， 北 斗 官方 正式 对 外 公布 北斗 坐标 系 BDCS (2019v01)， 它 是 通过 2019 年 
1 月 1 日 至 3 月 31 日 ,全球 均匀 分 布 120 多 个 台 站 (包括 IGS 台 站 和 北斗 区 域 监测 站 ) 的 
GPS 观测 数据 实现 并 对 齐 ITRF2014， 精 度 为 毫米 级 ””。 由 于 目前 采用 的 观测 数据 时 间 
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较 短 (< 2.5 a)， 无 法 估计 精确 的 台 站 速度 ”， 因 此 BDCS 采用 高 频率 更 新 方式 (更 新 周期 
为 1 a)， 在 更 新 周期 内 对 各 个 监测 站 持续 进行 监测 ， 以 此 来 维持 BDCS 毫米 级 精度 。 
目前 无 论 是 CGCS2000 还 是 BDCS， 都 是 通过 GPS 观测 数据 实现 和 维持 ， 而 不 是 依托 
北斗 单 系统 观测 数据 ， 基 于 此 ， 本 文 首先 利用 2019 年 1 月 至 2020 年 4 月 北斗 区 域 监测 台 
站 、 中 国 大 陆 环境 构造 监测 网 络 ( 陆 态 网 ) 台 站 、 全 球 分 布 IGS 和 iGMAS 台 站 的 单 GPS 观 
测 数据 ， 实 现 全 球 参考 框架 并 对 齐 ITRF2014， 利 用 IGS 提供 的 单 天 台 站 坐标 对 计算 方法 、 
分 析 模 型 和 结果 的 可 靠 性 进行 验证 。 其 次 ， 采 用 相同 的 计算 方法 和 分 析 模 型 ， 利 用 2019 年 
10 月 至 2020 年 3 月 北斗 区 域 监测 站 、 全 球 均 匀 分 布 IGS 台 站 和 iGMAS 台 站 的 单 北 斗 系 统 
观测 数据 ， 初 步 尝试 获得 历 元 台 站 坐标 精确 估计 ， 实 现 BDCS 并 与 ITRF2014 对 齐 ， 利 用 
单 GPS 数据 解 算 的 台 站 坐标 和 IGS 提供 的 台 站 坐标 进行 精度 评估 。 


2 数据 处 理 方法 


本 文 利用 GAMIT/GLOBK 软件 名 进行 GNSS 数据 处 理 ， 根 据 台 站 数量 分 成 几 个 子 网 ， 
选取 一 定数 量 的 稳定 台 站 作为 子 网 中 的 公共 站 ， 每 个 子 网 对 卫星 精密 轨道 、 大 气 参数 、 极 
移 、UT1 和 坐标 等 参数 进行 估计 ， 获 得 台 站 坐标 的 单 天 松弛 解 。 对 于 IGS 台 站 ， 如 果 有 
ITRF2014 坐标 ， 初 坐标 采用 ITRF2014 坐标 ， 其 他 站 采用 单 点 定位 坐标 为 初 值 。IGS 核心 
站 的 先 验 精 度 为 5 cm， 其 他 站 先 验 精度 为 100 m。 所 有 子 网 单 天 松弛 解 算 之 后 ， 采 用 基于 
卡尔 曼 滤 波 的 瞬时 历 元 对 齐 ITRF 方法 ”， 对 每 个 历 元 所 有 子 网 的 单 天 松弛 解 通过 最 小 约 
束 方法 进行 综合 解 算 ， 获 得 台 站 历 元 坐标 精确 估计 ， 实 现 瞬 时 历 元 北斗 坐标 系 ， 通 过 约束 
IGS 核心 站 的 ITRF2014 坐标 ， 将 台 站 坐标 对 齐 到 ITRF2014。 最 小 约束 方法 即 让 参考 台 站 


TGS 全 让 iGMAS 台 站 | 坐标 相对 先 验 坐 标的 坐标 差 之 间 在 最 小 二 乘法 

AOE) CUMAE L AEE] CRUSHES. REALE. P 
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| 单 天 全 球 网 松 约束 解 算 | 个 同 参考 框架 之 间 的 转换 通常 采用 Helmert 

es] 变换 方法 ， 对 于 两 个 不 同 参考 框架 之 间 的 转换 


采用 7/14 参数 还 是 6/12 参数 ， 目 前 国际 上 仍 


对 齐 ITRF v [34,35 、 7 

| ERGO, BI MRR 7 参 
HEF 数 ， 为 了 实现 各 系统 之 间 的 兼容 互 操作 ， 本 文 
+ 采用 7 SB. WETER ITRF 流程 如 


图 1 所 示 ，GAMIT/GLOBK 中 采用 的 分 析 模 
型 如 表 1 所 示 。 


图 1 数据 处 理 流程 


3 GPS 数据 实现 全 球 参 考 框 架 对 齐 ITRF2014 的 精度 评估 


为 了 对 第 4 章 采 用 单 BDS 数据 计算 结果 的 准确 性 和 可 靠 性 进行 验证 ， 本 章 首 先 利用 全 


球 均匀 分 布 台 站 的 单 GPS 数据 估算 历 元 台 站 精确 坐标 和 人 台 站 速度 ， 实 现 全 球 参考 框架 对 齐 
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表 1 数据 处 理 采 用 的 模型 和 估计 参数 

模型 参数 

观测 值 GNSS 卫星 双 差 无 电离 层 组 合 观测 值 (LC) BDS: B1iB3i; GPS: L1L2 

高 度 截止 角 10° 

采样 率 30s 

卫星 天 线 相位 中 心 ”IGS14.atx 

对 流 层 IERS2010 convention 

高 阶 电 离 层 二 阶 电离 层 

固体 消 IERS2010 Con 

AEN FES2004 model 

极 漳 IERS2010 convenon a 

相对 论 IERS 2010 convention | 

NN 体 JPL DE403 

重力 场 EGM2008 10x10 

太阳 辐射 压 ECOME 
估计 参数 

卫星 轨道 位 置 和 速度 ， 每 颗 星 每 个 弧 段 一 组 

太阳 光 压 每 颗 星 每 个 弧 段 一 组 

大 气 时 延 2h 一 组 

大 气 梯度 南北 、 东 西方 向 每 12 h 一 组 

台 站 坐标 每 天 


> ITRF2014， 并 通过 IGS 提供 台 站 坐标 序列 进行 精度 评估 。 
3.1 ”数据 及 台 站 分 布 
采用 2019 年 1 月 1 日 起 至 2020 年 4 月 20 日 220 个 台 站 的 GPS 数据， 包括 北斗 区 域 
监测 站 、 全 球 均匀 分 布 的 IGS 站 、 中 国 大 陆 环境 构造 监测 网 络 (CMONOC 陆 态 网 ) 台 站 和 
iGMAS 台 站 ， 台 站 分 布 如 图 2 所 示 。 为 了 对 区 域 监测 站 进行 长 期 监测 ， 在 每 个 区 域 监测 站 
ia 附近 选取 陆 态 网 和 iGMAS 并 址 台 站 同时 参与 解 算 。 在 BDCS (2019V01) 版 本 中 采用 的 观测 
台 站 数量 基础 上 增加 了 能 接受 北斗 三 号 卫星 观测 数据 的 IGS 台 站 ， 精 密 轨道 和 钟 差 文 件 来 
HF IGS (https://cddis.nasa.gov/). 48 220 个 台 站 分 成 6 个 子 网 ， 每 个 子 网 中 公共 台 站 有 
40 “SIGS 稳定 核心 站 ， 采 用 的 松 约束 解 算 方法 、 模 型 、 约 束 标准 以 及 对 齐 ITRF2014 方法 
如 第 2 章 介 绍 。 图 3 给 出 每 天 有 效 的 台 站 数量 、 对 齐 ITRF2014 的 约束 台 站 数量 以 及 观测 
卫星 数量 ， 可 以 看 出 每 天 观测 台 站 数量 在 2019 年 初 为 160 个 左右 ， 到 2020 年 3 月 增加 到 
200 个 台 站 左右 ， 从 2019 年 8 月 开始 ， 台 站 数量 增加 。 这 是 因为 在 2019 年 初期 部 分 北斗 区 
域 台 站 运行 尚 不 稳定 ， 随 着 时 间 累 积 ， 区 域 台 站 运行 稳定 后 ， 参 与 解 算 的 区 域 台 站 及 观测 数 
据 增 多 ， 部 分 历 元 有 效 观测 台 站 数量 少 于 150 个 也 是 由 区 域 台 站 数量 少 所 引起 的 。 每 天 有 

效 对 齐 ITRF2014 的 约束 台 站 数量 为 25 ~ 40 个 ， 观 测 卫 星 数 量 为 30, 31 或 者 32 个 。 
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图 2 全球 IGS 台 站 (紫色 )、iGMAS 台 站 (HE) 和 陆 态 网 台 站 (红色 ) 分 布 图 
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图 3 每 天 有 效 的 台 站 数量 (上 图 黑色 )、 对 齐 ITRF2014 的 约束 台 站 数量 (上 图 红色 ) 和 有 效 的 卫星 数量 


3.2 “” 台 站 坐标 精度 评估 
3.2.1 并 址 站 坐标 时 间 序 列 

从 260 个 陆 态 网 GNSS 观测 台 站 中 选取 与 区 域 监测 站 距离 最 近 的 台 站 作为 并 址 站 ， 通 
过 并 址 台 站 坐标 时 间 序 列 可 以 对 区 域 站 的 长 期 稳定 性 进行 评估 。 以 北京 地 区 为 例 ， 该 地 区 有 
4 个 并 址 台 站 ， 图 4 为 各 台 站 坐标 时 间 序列 ， 图 中 BJFS (SHAO) 为 本 次 解 算 BJFS 台 站 ， 
BJSH 为 陆 态 网 台 站 ，BJF1 NiGMAS 台 站 ，3012 为 区 域 监测 站 ，BJFS (IGS) 为 IGS 解 
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算 的 BJFS 站 。 为 方便 比较 ， 图 中 每 个 台 站 水 平方 向 偏 移 10 mm, FETT Al ts 20 mm。 
由 图 可 知 ，5 个 台 站 坐标 时 间 序 列 水 平 与 高 程 方向 具有 较 好 一 致 性 ， 水 平方 向 存在 明显 的 趋 
势 性 变化 。 图 5 为 内 蒙古 地 区 4 个 并 址 台 站 坐标 时 间 序 列 ， 高 程 方向 共有 非 线性 特性 。 


*BJFS(SHAO) 
= .BJSH 
ra -BJF1 
a -3012 

*BJFS(IGS) 
i 
g 
E 
局 

50 L 4. i 4 4 上 
2019 2019.2 2019.4 2019.6 2019.8 2020 2020.2 2020.4 
年 份 
注 : E, N, U 分 别 代表 台 站 坐标 东西 、 南 北 、 高 程 的 量 。 
图 4 北京 地 区 并 址 台 站 坐标 时 间 序 列 
-3361 
g + 3362 
Ë +3363 
S -HLHG 
i 
g 
E 
S 


0 1 1 1 1 1 1 

2019 2019.2 2019.4 20196 2019.8 2020 2020.2 2020.4 
年 份 

5 ”内 蒙古 地 区 4 个 并 址 台 站 坐标 时 间 序 列 


3.2.2 ” 台 站 坐标 残 差 RMS 

采用 2019 年 1 月 1 日 至 2020 年 4 月 19 日 IGS 提供 的 天 解 坐标 值 (https://cddis.nasa. 
gov/) 对 本 次 解 算 的 98 SIGS 台 站 坐标 精度 进行 检验 。 图 6a) 为 解 算 的 ABPO 台 站 与 IGS 
提供 的 天 解 坐 标 时 间 序 列 图 ， 图 6b) 为 两 种 序列 坐标 差异 图 ， 两 种 台 站 坐标 时 间 序 列 吻合 
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性 较 好 。 图 7 为 98 个 台 站 坐标 在 E N 和 U 三 个 方向 差异 的 RMS Zit, sta REGU, 
RMS 均值 在 互 方 向 为 1.8 mm，N 方向 为 1.8 mm, U 方向 为 43 mm， 结 果 表 明 采 用 GPS 
数据 解 算 IGS 台 站 坐标 重复 率 优 于 1 cm. 


er = eee q ot ra 
ws 
2 J 
= ] aw? fe 
ESS EA P 。 


? | z 。 

2019.2 2019.6 2020 2020.4 i 2019.2 2019.6 2020 2020.4 
“eth Ef 
a) b) 


TE: a) 为 解 算 坐 标 ( 蓝 色 ) 与 IGS 提供 天 解 坐 标 (红色 ) 时 间 序列 ，b) 为 两 种 序列 坐标 差异 时 间 序 列 。 
图 6 ABPO 台 站 坐标 与 IGS 提供 台 站 坐标 时 间 序 列 


sta 


图 7 98 个 IGS 台 站 坐标 与 IGS 提供 天 解 坐标 差异 的 RMS 
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对 于 北斗 监测 台 站 ， 由 于 观测 数据 时 间 跨 度 短 (小 于 2.5 a)， 采 用 扣除 线性 趋势 项 方法 
对 台 站 坐标 精度 进行 分 析 。 图 8 为 区 域 监测 台 站 (3172) 的 坐标 时 间 序 列 。 图 9 为 35 个 区 域 
监测 站 扣除 趋势 项 后 残 差 的 RMS Kit, RMS 均值 已 方向 为 1.9 mm, N 方向 为 2.2 mm, 
U 方向 为 8.2 mm， 结 果 表明 采用 GPS 数据 解 算 区 域 监测 站 坐标 重复 率 优 于 1 cm. 
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N/mm 
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U/mm 


-20+ pi 
2019.2 2019.6 2020 2020.4 2019.2 2019.6 2020 2020.4 
年 份 年 从 
a) b) 

HE: a) 解 算 台 站 坐标 时 间 序 列 (红色 ) 和 拟 合 趋势 项 时 间 序 列 ( 蓝 色 ); b) 扣除 趋势 项 后 残 差 时 间 序 列 。 


图 8 采用 GPS 数据 解 算 区 域 台 站 坐标 时 间 序 列 


RMS,/mm RMS,/mm 


RMS,/mm 


0 5 10 15 20 25 30 35 


sta 


9 35 个 区 域 监测 站 坐标 时 间 序 列 扣除 趋势 项 后 残 差 的 RMS 
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3.3” 台 站 速度 场 

本 次 采用 的 GPS 观测 数据 时 间 跨 度 较 短 (470 d)， 因 此 只 考虑 线性 项 ， 通 过 线性 拟 合 得 
到 台 站 速度 。IGb14 台 站 速度 采用 20 多 年 GPS 观测 数据 进行 解 算 ， 并 考虑 周年 、 半 周年 、 
地 震 、 接 收 机 或 者 天 线 变化 等 各 种 因素 进行 非 线性 建 模 ， 最 终 得 到 精确 台 站 速度 ， 因 此 本 文 
以 IGb14 提供 的 台 站 速度 为 标准 对 解 算 的 98 个 IGS 台 站 速度 进行 精度 评估 。 图 10 为 98 个 
台 站 速度 差异 图 ，vsHAo 为 本 文 拟 合 的 台 站 速度 ，wrcs 为 对 IGS 提供 的 台 站 天 解 坐标 序列 
进行 线性 拟 合 得 到 的 台 站 速度 ，wcba 为 IGb14 提供 的 台 站 速度 。 从 图 中 可 以 看 出 ，vsHAo 
与 viGs 差异 的 RMS 在 2 mmy/a 以 内 ， 与 viGb14 差异 水 平方 向 在 2 mm/a WA, FEAT 
约 3 mm/a。 图 10b) F KOUR 台 站 (sta = 40) 和 OUS2 台 站 (sta = 60) 与 IGb14 速度 水 
平方 向 差异 较 大 ， 主 要 是 由 于 KOUR 在 2019 年 10 月 9 日 OUS2 在 2019 年 9 月 25 日 分 
别 更 换 接收 机 引起 。 
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10 98 个 IGS 台 站 速度 差异 


对 于 北斗 区 域 监测 站 ， 因 为 在 其 附近 选取 其 他 观测 网 络 的 并 址 台 站 的 同时 也 参与 解 算 ， 
所 以 通过 并 址 台 站 速度 对 区 域 监测 台 站 速度 进行 精度 评估 。 表 2 所 示 为 北京 地 区 各 并 址 台 站 
的 速度 值 ， 表 中 BJSH 为 解 算 的 陆 态 网 台 站 速度 ，BJF1 为 解 算 的 iGMAS 台 站 速度 ，BJFS 
(IGS) 为 对 IGS 提供 天 解 坐标 时 间 序 列 进行 线性 拟 合 得 到 的 BJFS 台 站 速度 ，3012 为 解 算 
的 北斗 区 域 监测 台 站 速度 ，BJFS (SHAO) 为 解 算 的 BJFS 台 站 速度 ，BJFS (ITRF2014) 为 


AT 
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ITRF2014 提供 的 BJFS 站 速度 ，BJFS (IGb14) 为 IGb14 提供 的 BJFS 台 站 速度 。 从 BJFS 
(SHAO) 分 别 与 BJFS (ITRF2014) 和 BJFS (IGb14) 速度 的 比较 可 以 看 出 ， 观 测 数据 时 间 
跨度 较 短 ， 对 速度 估计 的 精度 影响 水 平方 向 约 为 1 mm/a， 垂 直方 向 约 为 2 mm/a。BJFS 
(SHAO) 与 BJFS (IGS) 相 比 ，N 方向 速度 差异 约 1 mm/a, E 和 U 方向 差异 约 2 mm/a。 
3012 站 和 BIFL 站 由 于 部 分 观测 数据 缺失 以 及 未 对 台 站 序列 中 存在 的 非 线性 运动 进行 分 析 
建 模 ， 导 致 两 个 台 站 速度 在 包 和 U 方向 精度 相对 较 差 。 


表 2 北京 地 区 并 址 台 站 速度 


台 站 vw/mma sl! vp/mma! vu/mm-a ! 

BJSH —10.3 28.4 3.1 

BJF1 —11.0 26.5 一 13.8 

3012 一 12.0 25.1 一 8.3 
BJFS (SHAO) —11.8 29.3 4.4 
BJFS (IGS) —11.6 31.7 1.9 
BJFS (ITRF2014) —10.6 30.6 2.2 
BJFS (IGb14) —10.5 30.7 2.2 


图 11 为 中 国 区 域 并 址 台 站 的 水 平方 向 速度 场 图 ( 蓝 色 和 红色 分 别 代表 两 个 不 同 观测 网 
络 台 站 )， 从 图 中 可 以 看 出 ， 大 部 分 台 站 水 平方 向 速度 变化 趋势 一 致 性 较 好 ， 部 分 台 站 差异 
相对 较 大 是 由 观测 数据 时 间 跨 度 太 短 引 起 。 例 如 图 中 左上 方 蓝 色 框 内 的 台 站 速度 ， 红 第 头 代 
表 的 台 站 观测 数据 时 间 跨 度 有 450 d， 而 蓝 箭头 代表 的 台 站 观测 数据 时 间 跨 度 只 有 105 do 


72 84 96 108 120 132 
经 度 /(?) 


11 ”中 国 区 域 并 址 台 站 水 平方 向 速度 场 


从 台 站 速度 比较 可 以 看 出 ， 此 次 解 算 的 台 站 速度 场 与 IGb14 和 ITRF2014 台 站 速度 有 
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差异 ， 一 个 主要 原因 是 观测 数据 的 时 间 跨 度 短 (470 d)， 在 计算 台 站 速度 时 未 考虑 季节 项 、 
构造 运动 产生 的 跳跃 项 、 非 构造 运动 产生 的 跳跃 项 、 震 后 形变 和 其 他 形变 等 。Blewitt 和 
Lavelle” 研究 表明 ， 观 测 数据 时 间 长 度 至 少 2.5 a， 才 能 得 到 精确 的 台 站 速度 。 因 此 随 着 观 
测 数据 增加 ， 精 确 考 虑 每 个 台 站 的 各 种 非 线 性 因素 后 得 到 的 速度 精度 会 更 高 。 


4 北斗 单 系统 数据 实现 BDCS 对 齐 ITRF2014 的 精度 评估 


KI 


采用 2019 年 10 H 1 HÆ 2020 年 3 H 19 日 全 球 分 布 台 站 接收 到 的 北斗 观测 数据 实现 


BDCS 并 对 齐 ITRF2014， 并 与 第 3 BH GPS 解 算 结果 及 IGS 提供 的 台 站 坐标 进行 比较 


分 析 。 


4.1 数据 及 台 站 分 布 
采用 北斗 单 系统 观测 数据 进行 解 算 时 ， 轨 道 和 钟 差 产品 来 自 于 iGMAS (http://www. 
igmas.org/)。 由 于 在 本 文 计算 时 ，2019 年 10 月 1 日 之 前 没有 北斗 BII 和 了 B3I 信号 的 精密 钟 


差 产品 ， 


因此 选取 从 2019 年 10 月 1 日 至 2020 年 3 月 19 日 能 够 接收 北斗 B1I 和 B3I 频率 


观测 数据 的 全 球 209 个 台 站 ， 包 括 北 斗 区 域 监测 站 、 全 球 均匀 分 布 IGS 核心 台 站 和 iGMAS 
全 球 分 布 台 站 ， 如 图 12 所 示 。 


12 IGS 台 站 (紫色 ) 和 iGMAS £% (HE) 分 布 图 


把 209 个 台 站 分 成 5 个 子 网 ， 从 IGS14 核心 站 中 选取 能 够 接收 北斗 三 号 观测 数据 的 40 
个 台 站 作为 每 个 子 网 中 的 公共 台 站 ， 每 个 子 网 的 松 约束 解 算 方法 、 采 用 的 模型 、 约 束 标准 以 
及 对 齐 ITRF2014 方法 如 第 2 章 介 绍 。 每 个 历 元 有 效 台 站 数量 、 对 齐 ITRF2014 WARE 
站 数量 以 及 观测 卫星 数量 如 图 13 所 示 ，2019 年 10 月 每 个 历 元 有 效 观 测 台 站 数量 约 150 个 ， 


到 2020 


年 3 月 增加 到 200 个 左右 。 对 齐 ITRF2014 的 约束 台 站 数量 从 25 个 左右 增加 至 40 
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个 左右 ， 观 测 卫星 数量 从 20 个 左右 增加 至 30 个 左右 。 统 计 表 明 ， 随 着 时 间 积 累 ， 全 球 能 
接收 北斗 三 号 卫星 观测 数据 的 台 站 越 来 越 多 ， 可 用 的 有 效 观测 数据 和 卫星 数量 也 越 来 越 多 。 


200 T en pa, 30 T | r T ' ret 
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a) b) 


VE: a) 黑色 表示 每 个 历 元 有 效 台 站 数量 ， 红 色 表 示 对 齐 ITRF2014 的 约束 台 站 数量 ，b) 卫星 数量 。 


图 13 “” 台 站 数量 和 卫星 数量 随时 间 分 布 


4.2 ” 台 站 坐标 精度 评估 

4.2.1 区 域 监测 站 坐标 残 差 RMS 

通过 台 站 坐标 时 间 序列 消除 线性 项 后 残 差 的 RMS 对 31 个 区 域 监测 站 坐标 精度 进行 分 
Ff。 图 14a) 为 台 站 (3172) 坐标 时 间 序 列 以 及 拟 合 的 线性 趋势 项 ， 图 b) 为 消除 趋势 项 后 的 
标 残 差 时 间 序 列 图 。 图 15 为 31 个 区 域 监测 站 消除 线性 项 后 残 差 的 RMS 统计 图 ， 可 以 
看 出 大 部 分 台 站 水 平方 向 RMS 小 于 5 mm， 垂 直方 向 RMS 小 于 8 mm，RMS WE 方 
HX 4.4mm, N JHA 3.4mm, U 方向 为 7.9 mm， 结 果 表 明 北 斗 单 系统 数据 解 算 的 区 域 
监测 站 坐标 重复 率 优 于 1 cm 
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14 ”监测 站 (3172) 坐标 时 间 序 列 以 及 消除 线性 项 后 残 差 图 
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图 15 北斗 单 系统 数据 解 算 31 个 区 域 站 坐标 时 间 序 列 消除 趋势 项 后 残 差 的 RMS 


4.2.2 IGS 台 站 坐标 残 差 RMS 

利用 IGS 提供 的 台 站 坐标 对 解 算 的 IGS 台 站 坐标 精度 进行 评估 。 图 16 列 出 ABPO & 
站 坐标 时 间 序 列 、IGS 提供 天 解 台 站 坐标 时 间 序列 以 及 两 种 台 站 坐标 差异 的 时 间 序 列 图 ，a) 
中 两 种 时 间 序 列 趋势 性 一 致 ， 但 是 采用 BDS 数据 解 算 的 台 站 坐标 时 间 序 列 基 弥散 度 明显 全 
大 。 图 17 为 解 算 的 126 个 IGS 台 站 坐标 与 IGS 提供 台 站 坐标 差异 RMS, EFA RMS 均 
值 5.96 mm, N 方向 RMS 均值 4.5 mm， 重 直方 向 RMS 均值 为 12.5 mm， 结 果 表 明 北 斗 


单 系统 观测 数据 解 算 IGS 台 站 坐标 重复 率 约 为 1 cm. 
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注 : a) BDS 数据 解 算 台 站 ( 蓝 色 ) SIGS 天 解 台 站 (红色 ) 的 坐标 时 间 序 列 ; b) 两 种 坐标 差异 时 间 序 列 。 
图 16 ABPO 台 站 坐标 时 间 序 列 图 
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17 BDS 数据 解 算 126 SIGS 台 站 坐标 与 IGS 天 解 坐 标 差 异 的 RMS 
除 采 用 IGS 提供 的 台 站 坐标 对 解 算 的 IGS 台 站 坐标 进行 精度 评估 外 ， 我 们 用 第 3 章 


中 单 GPS 数据 解 算 的 台 站 坐标 对 BDS 数据 解 算 结果 进行 精度 评估 。 图 18 A BDS 数据 和 
GPS 数据 分 别 解 算 ABPO 台 站 坐标 时 间 序 列 及 两 者 差异 图 ，a) 中 两 种 时 间 序 列 趋势 性 一 


致 ， 但 是 BDS 数据 解 算 的 台 站 坐标 时 间 序 列 弥散 度 偏 大 。 图 19 为 统计 BDS 数据 与 GPS 


数据 分 别 解 算 的 106 个 台 站 坐标 差异 (IGS 台 站 和 区 域 监测 台 站 )， 大 部 分 台 站 坐标 差异 的 
RMS 在 水 平和 垂直 方向 分 别 小 于 1 cm 和 2 cm, RMS WEE WAN 7.0mm, N 方向 
A48mm, U 方向 为 13.8 mm， 结 果 表 明 单 BDS 与 单 GPS 数据 解 算 台 站 坐标 差异 优 于 
1.5 cm. 
20 20 
g ~ te g : oe 
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a) b) 
iE: a) BDS 解 算 ( 蓝 色 ) 与 GPS 解 算 (红色 ) 时 间 序列 ，b) 两 者 差异 时 间 序列 。 
18 #BDS 数据 与 单 GPS 数据 解 算 ABPO 台 站 坐标 时 间 序 列 图 
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图 19 SA BDS 数据 与 单 GPS 解 算 台 站 坐标 差异 的 RMS 


如 前 所 述 ， 图 16 和 18 显示 BDS 数据 解 算 的 台 站 坐标 时 间 序 列 无 论 与 GPS 数据 解 算 
坐标 时 间 序 列 还 是 IGS 提供 台 站 天 解 坐标 时 间 序 列 比较 ， 趋 势 性 吻合 得 都 较 好 ， 但 是 BDS 
解 算 台 站 坐标 时 间 序 列 弥散 度 偏 大 。 为 了 解释 这 一 原因 ， 采 用 几 组 不 同 观测 数据 和 精密 轨 
道 产 品 解 算 台 站 坐标 进行 比较 分 析 ， 如 表 3 所 示 。casel 采用 209 个 台 站 的 BDS 观测 数据 ， 
BDS 精密 轨道 产品 来 自 于 iGMAS， 即 本 节 解 算 方 法 ，case2 采用 与 case 1 相同 的 台 站 和 相 
同时 间 跨 度 的 GPS 观测 数据 ，GPS 精密 轨道 产品 来 自 于 iGMAS; case3 采用 与 casel 相同 
的 台 站 和 相同 时 间 跨 度 的 GPS 观测 数据 ，GPS 精密 轨道 产品 来 自 于 IGS; cased 为 IGS 提 
供 单 天 台 站 坐标 解 (SINEX 格式 )， 作 为 参考 标准 对 其 他 3 组 解 算 结果 进行 评估 。 


表 3 4 组 不 同 试验 介绍 
case “ 台 站 数量 观测 数据 时 间 跨 度 观测 数据 类 型 ”精密 轨道 产品 


1 209 2019/10/01 一 2020/03/19 BDS iGMAS 
2 209 2019/10/01 一 2020/03/19 GPS iGMAS 
3 209 2019/10/01 一 2020/03/19 GPS IGS 
4 ”大 于 300 2019/10/01 一 2020/03/19 GPS IGS 


图 20 为 case2, case3 和 case4 分 别 解 算 的 ABPO 台 站 坐标 时 间 序 列 ， 图 中 case2 和 
case3 的 台 站 坐标 时 间 序 列 弥散 度 都 略 大 于 case4。 图 21a) 为 case2 和 case3 解 算 的 140 个 
IGS 台 站 坐标 时 间 序 列 差 异 的 RMS 统计 图 ，RAMS 均值 在 水 平方 向 小 于 4 mm， 垂 直方 
向 为 5.3 mm。case2 和 case3 采用 相同 的 GPS 观测 数据 ， 精 密 轨道 分 别 来 自 于 iGMAS 和 
IGS， 结 果 表 明 iGMAS 的 轨道 产品 与 IGS 轨道 产品 精度 相当 ， 这 与 iGMAS 评估 结果 一 
gc", R 21b) 为 case3 和 cased 解 算 的 125 个 IGS 台 站 坐标 时 间 序 列 差异 的 RMS 统计 
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Al, RMS 均值 水 平方 向 小 于 4mm， 重 直方 向 为 6.2 mm。case3 和 cased 采用 相同 GPS 观 
测 数据 和 IGS 精密 轨道 ，IGS 选取 的 台 站 相对 更 多 ， 各 分 析 中 心 数据 处 理 时 同时 解 算 卫 星 
轨道 和 钟 差 ， 通 过 IGS 分 析 中 心 综合 剔除 粗 差 ， 因 此 cased 台 站 坐标 时 间 序 列 弥散 度 会 相 
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20 case2 (红色 )、case3 (HE) 和 case4 (RE) 分 别 解 算 的 ABPO 台 站 坐标 时 间 序 列 
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21 case2 和 case3 解 算 的 140 SIGS 台 站 坐标 时 间 序 列 互 差 的 RMS 统计 


图 22a) 为 casel 和 case2 分 别 解 算 的 ABPO 


台 站 坐标 时 间 序 列 ， 两 种 时 间 序 列 弥 散 度 


也 同样 大 于 cased 的 台 站 坐标 时 间 序 列 ，casel 与 cased 的 弥散 度 差 异 相 对 更 大 。casel 和 
case2 分 别 采用 BDS 和 GPS 观测 数据 ， 精 密 轨 道 都 来 自 iGMAS， 因 此 两 者 差别 主要 来 自 
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不 同 的 卫星 系统 数据 和 相应 的 分 析 模 型 。 图 22b) 为 casel 和 case2 解 算 的 140 个 IGS 台 站 
坐标 差异 的 RMS Site, RMS 均值 水 平方 向 7 mm 以 内 ， 垂 直方 向 为 13.7 mm, HAE 
斗 系统 分 析 模 型 、 观 测 数据 质量 、 卫 星 履 盖 率 以 及 可 用 的 卫星 数 相 对 较 少 等 是 引起 台 站 坐标 
时 间 序 列 弥 散 度 大 的 另 一 个 原因 。 随 着 以 后 北斗 观测 数据 以 及 可 用 卫星 数量 增加 ， 我 们 将 对 


具体 原因 进行 更 详细 深入 的 分 析 。 


> 


RMS,/mm RMS,/mm RMS,/mm 


2019.8 2019.9 2020 2020.1 2020.2 
年 份 
a) 


22 a) casel (21), case2 (HE) 和 cased (绿色 ) RA ABPO 台 站 坐标 时 间 序 列 ; b) casel 和 case2 
解 算 IGS 台 站 坐标 差异 的 RMS 


4.2.3 赫 尔 默 特 转换 参数 

对 某 一 参考 框架 与 ITRF 之 间 的 赫 尔 默 特 转换 参数 进行 分 析 , 可 以 探究 参考 框架 对 齐 
ITRF 的 自 洽 性 和 稳定 性 。 图 23 为 采用 BDS 数据 和 GPS 数据 ， 通 过 约束 IGS 核心 台 站 ， 
分 别 得 到 的 每 个 历 元 参考 框架 对 齐 ITRF2014 的 7 参数 时 间 序 列 。 图 23a) H Te, Ty M T, 
分 别 为 平移 参数 ， 最 下 面 图 为 尺度 时 间 序 列 ， 图 23b) 为 3 个 旋转 参数 时 间 序 列 图 。 从 图 中 
看 出 单 GPS 数据 解 算 的 7 参数 时 间 序 列 更 加 平稳 ，BDS 数据 解 算 的 时 间 序 列 弥 散 度 偏 大 ， 
尤其 在 2019 年 10 月 和 11 月 期 间 ， 这 可 能 是 由 于 北斗 三 号 系统 运行 初期 ， 可 用 卫星 数量 、 
观测 数据 相对 较 少 等 原因 引起 ， 也 可 能 是 由 于 目前 北斗 导航 系统 中 可 能 存在 的 误差 ， 例 如 卫 
星相 位 中 心 模型 、 卫 星 轨道 光 压 模型 等 引起 的 不 稳定 误差 进入 7 参数 中 。 随 着 北斗 系统 更 加 
稳定 运行 ， 积 累 的 观测 数据 更 长 ， 我 们 将 对 7 个 转换 参数 进行 更 加 详细 精确 的 分 析 。 


本 文采 用 2019 年 10 月 1 日 至 2020 年 3 月 19 日 北斗 区 域 监测 站 和 IGS 台 站 的 北斗 
观测 数据 ， 初 步 尝 试 实现 历 元 BDCS 对 齐 ITRF2014， 估 计 每 个 历 元 台 站 精确 坐标 。 利 用 
2019 年 1 月 1 日 至 2020 Æ 4 H 20 H GPS 观测 数据 的 计算 结果 以 及 IGS 提供 的 天 解 台 站 
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图 23 RA BDS 数据 ( 蓝 色 ) 和 GPS 数据 (红色 ) 分 别 实现 每 个 历 元 参考 框架 对 齐 ITRF2014 的 7 参数 


时 间 序 列 


这 


坐标 时 间 序 列 进行 台 站 坐标 精度 评估 ， 上 所 得 


主要 结果 如 下 : 


(1) 单 GPS 数据 解 算 IGS 台 站 和 区 域 监测 站 坐标 重复 率 优 于 1 cm 


(2) 北斗 单 系统 数据 解 算 IGS 


台 站 坐标 与 IGS 天 解 台 站 坐标 差异 的 RMS 


均值 小 于 


1cm， 与 GPS 结果 互 差 RMS 均值 小 于 1.5 cm， 台 站 坐标 时 间 序 列 变化 趋势 相 一 致 ， 北 
斗 单 系统 数据 解 算 的 区 域 监测 站 坐标 时 间 序 列 ， 扣 除 线性 速度 之 后 残 差 的 RMS 


1 cmo 


(3) 北斗 单 系统 数据 解 算 台 站 坐标 时 间 序 列 药 
BDS 卫星 数据 质量 等 原因 有 关 。 北斗 身 
时 间 序 列 弥散 度 偏 大 ， 可 能 与 北斗 导航 系统 中 存在 的 误差 和 可 用 卫星 数量 


等 原因 有 关 。 


(4) 利用 GPS 数据 初步 实现 区 域 台 


稳定 的 参考 框架 
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均值 优 于 


弥散 度 偏 大 ， 这 可 能 与 框架 台 站 分 布 以 及 
和 系统 数据 实现 BDCS 对 齐 ITRF2014 的 历 元 7 参数 
及 观测 数据 量 少 


站 初步 速度 场 解 算 ， 随 着 数据 的 累积 ， 有 望 利 用 北 
斗 单 系统 数据 实现 台 站 的 速度 高 精度 估计 。 


(5) 随 着 BDS 的 稳定 运行 ， 后 续 有 望 综合 利用 BDS, GPS 和 Galileo 获得 精度 更 高 、 更 


i GPS P iGMAS 提供 精密 产品 、BDS 
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Accuracy Evaluation of Aligning BeiDou Coordinate System 
to ITRF2014 Based on BeiDou Observations 


QU Wei-jing!, LIU Li?, ZHOU Shan-shi!, XU Jun-yi?, YANG Yu-fei?, 
HU Xiao-gong!, GUO Rui?, WU Bin', DONG Ming’, 
TAND Cheng-pan', LI Xiao-jie?, LIU Shu-ai? 
(1. Shanghai Astronomical Observatory, Chinese Academy of Science, Shanghai 200030, China; 2. 


Beijing Satellite Navigation Center, Beijing 100092, China; 3. Beijing Institute of Tracking and T- 
elecommunication Technology, Beijing 100094, China) 


Abstract: Based on GPS observations, the BeiDou coordinate system (BDCS) (2019v01) 
aligned to ITRF2014 was released in December 2019 and the alignment accuracy is at mil- 
limeter level. With the continuous accumulation of BeiDou Satellite Navigation System 
(BDS) observations, it should be used to the alignment and maintenance of BDCS. In this 
paper, we align BDCS to ITRF2014 using BDS observations from October 1, 2019 to March 
19, 2020. The daily station coordinates obtained by using GPS observations from global 
distribution stations and that from the International GPS Service (IGS) from January 1, 
2019 to April 19, 2020 are used to evaluate the accuracy of the station coordinates estimated 
by using BDS observations. The average root-mean-square (RMS) of the daily IGS station 
coordinate differences between using BDS observations and from IGS analysis center is less 
than 1 cm. The average RMS of the daily station coordinate differences between using G- 
PS and BDS observations is less than 1.5 cm. The scatter of station coordinates solutions 
using BDS observation is relatively large, which is preliminary explained in this paper. The 
velocities of the regional monitoring stations are estimated using GPS observations of 470 
days. With the stable operation of BDS and continued accumulation of BDS observations, 
we expect to obtain higher accuracy of station coordinates, velocity, and the alignment to 
ITRF. 


Key words: BeiDou Satellite Navigation System; BeiDou coordinate system; International 


celestial reference frame 


